Different receptor models show differences in ligand binding strength and location: a computational drug screening for the tick-borne encephalitis virus.

The tick-borne encephalitis virus (TBE) is a neurotrophic disease that has spread more rapidly throughout Europe and Asia in the past few years. At the same time, no cure or specific therapy is known to battle the illness apart from vaccination. To find a pharmacologically relevant drug, a computer-aided drug screening was initiated. Such a procedure probes a possible binding of a drug to the RNA Polymerase of TBE. The crystal structure of the receptor, however, includes missing and partially modeled regions, which rendered the structure incomplete and of questionable use for a thorough drug screening procedure. The quality of the receptor model was addressed by studying three putative structures created. We show that the choice of receptor models greatly influences the binding affinity of potential drug molecules and that the binding location could also be significantly impacted. We demonstrate that some drug candidates are unsuitable for one model but show decent results for another. Without any prejudice on the three employed receptor models, the study reveals the imperative need to investigate the receptor structure before drug binding is probed whether experimentally or computationally.

Structural Rearrangements of Pigeon Cryptochrome 4 Undergoing a Complete Redox Cycle.

Cryptochrome is currently the major contender of a protein to underpin magnetoreception, the ability to sense the Earth's magnetic field. Among various types of cryptochromes, cryptochrome 4 has been identified as the likely magnetoreceptor in migratory birds. All-atom molecular dynamics (MD) studies have offered first insights into the structural dynamics of cryptochrome but are limited to a short time scale due to large computational demands. Here, we employ coarse-grained MD simulations to investigate the emergence of long-lived states and conformational changes in pigeon cryptochrome 4. Our coarse-grained simulations complete the picture by permitting observation on a significantly longer time scale. We observe conformational transitions in the phosphate-binding loop of pigeon cryptochrome 4 upon activation and identify prominent motions in residues 440-460, suggesting a possible role as a signaling state of the protein or as a gated interaction site for forming protein complexes that might facilitate downstream processes. The findings highlight the importance of considering longer time scales in studying cryptochrome dynamics and magnetoreception. Coarse-grained MD simulations offer a valuable tool to unravel the complex behavior of cryptochrome proteins and shed new light on the mechanisms underlying their role in magnetoreception. Further exploration of these conformational changes and their functional implications may contribute to a deeper understanding of the molecular mechanisms of magnetoreception in birds.

Introducing the Automated Ligand Searcher.

The Automated Ligand Searcher (ALISE) is designed as an automated computational drug discovery tool. To approximate the binding free energy of ligands to a receptor, ALISE includes a three-stage workflow, with each stage involving an increasingly sophisticated computational method: molecular docking, molecular dynamics, and free energy perturbation, respectively. To narrow the number of potential ligands, poorly performing ligands are gradually segregated out. The performance and usability of ALISE are benchmarked for a case study containing known active ligands and decoys for the HIV protease. The example illustrates that ALISE filters the decoys successfully and demonstrates that the automation, comprehensiveness, and user-friendliness of the software make it a valuable tool for improved and faster drug development workflows.

Across atoms to crossing continents: Application of similarity measures to biological location data.

Biological processes involve movements across all measurable scales. Similarity measures can be applied to compare and analyze these movements but differ in how differences in movement are aggregated across space and time. The present study reviews frequently-used similarity measures, such as the Hausdorff distance, Fréchet distance, Dynamic Time Warping, and Longest Common Subsequence, jointly with several measures less used in biological applications (Wasserstein distance, weak Fréchet distance, and Kullback-Leibler divergence), and provides computational tools for each of them that may be used in computational biology. We illustrate the use of the selected similarity measures in diagnosing differences within two extremely contrasting sets of biological data, which, remarkably, may both be relevant for magnetic field perception by migratory birds. Specifically, we assess and discuss cryptochrome protein conformational dynamics and extreme migratory trajectories of songbirds between Alaska and Africa. We highlight how similarity measures contrast regarding computational complexity and discuss those which can be useful in noise elimination or, conversely, are sensitive to spatiotemporal scales.

Effects of Dynamical Degrees of Freedom on Magnetic Compass Sensitivity: A Comparison of Plant and Avian Cryptochromes.

The magnetic compass of migratory birds is thought to rely on a radical pair reaction inside the blue-light photoreceptor protein cryptochrome. The sensitivity of such a sensor to weak external magnetic fields is determined by a variety of magnetic interactions, including electron-nuclear hyperfine interactions. Here, we investigate the implications of thermal motion, focusing on fluctuations in the dihedral and librational angles of flavin adenine dinucleotide (FAD) and tryptophan (Trp) radicals in cryptochrome 4a from European robin (Erithacus rubecula, ErCry4a) and pigeon (Columba livia, ClCry4a) and cryptochrome 1 from the plant Arabidopsis thaliana (AtCry1). Molecular dynamics simulations and density functional theory-derived hyperfine interactions are used to calculate the quantum yield of radical pair recombination dependent on the direction of the geomagnetic field. This quantity and various dynamical parameters are compared for [FAD•- Trp•+] in ErCry4a, ClCry4a, and AtCry1, with TrpC or TrpD being the third and fourth components of the tryptophan triad/tetrad in the respective proteins. We find that (i) differences in the average dihedral angles in the radical pairs are small, (ii) the librational motions of TrpC•+ in the avian cryptochromes are appreciably smaller than in AtCry1, (iii) the rapid vibrational motions of the radicals leading to strong fluctuations in the hyperfine couplings affect the spin dynamics depending on the usage of instantaneous or time-averaged interactions. Future investigations of radical pair compass sensitivity should therefore not be based on single snapshots of the protein structure but should include the ensemble properties of the hyperfine interactions.

Computational Reconstruction and Analysis of Structural Models of Avian Cryptochrome 4.

A recent study by Xu et al. (Nature, 2021, 594, 535-540) provided strong evidence that cryptochrome 4 (Cry4) is a key protein to endow migratory birds with the magnetic compass sense. The investigation compared the magnetic field response of Cry4 from migratory and nonmigratory bird species and suggested that a difference in magnetic sensitivity could exist. This finding prompted an in-depth investigation into Cry4 protein differences on the structural and dynamic levels. In the present study, the pigeon Cry4 (ClCry4) crystal structure was used to reconstruct the missing avian Cry4 protein structures via homology modeling for carefully selected bird species. The reconstructed Cry4 structure from European robin, Eurasian blackcap, zebra finch, chicken, and pigeon were subsequently simulated dynamically and analyzed. The studied avian Cry4 structures show flexibility in analogous regions pointing to similar activation mechanisms and/or signaling interaction partners. It can be concluded that the experimentally recorded difference in the magnetic field sensitivity of Cry4 from different birds is unlikely to be due to solely intrinsic dynamics of the proteins but requires additional factors that have not yet been identified.

The Transition of Photoreceptor Guanylate Cyclase Type 1 to the Active State.

Membrane-bound guanylate cyclases (GCs), which synthesize the second messenger guanosine-3', 5'-cyclic monophosphate, differ in their activation modes to reach the active state. Hormone peptides bind to the extracellular domain in hormone-receptor-type GCs and trigger a conformational change in the intracellular, cytoplasmic part of the enzyme. Sensory GCs that are present in rod and cone photoreceptor cells have intracellular binding sites for regulatory Ca2+-sensor proteins, named guanylate-cyclase-activating proteins. A rotation model of activation involving an α-helix rotation was described as a common activation motif among hormone-receptor GCs. We tested whether the photoreceptor GC-E underwent an α-helix rotation when reaching the active state. We experimentally simulated such a transitory switch by integrating alanine residues close to the transmembrane region, and compared the effects of alanine integration with the point mutation V902L in GC-E. The V902L mutation is found in patients suffering from retinal cone-rod dystrophies, and leads to a constitutively active state of GC-E. We analyzed the enzymatic catalytic parameters of wild-type and mutant GC-E. Our data showed no involvement of an α-helix rotation when reaching the active state, indicating a difference in hormone receptor GCs. To characterize the protein conformations that represent the transition to the active state, we investigated the protein dynamics by using a computational approach based on all-atom molecular dynamics simulations. We detected a swinging movement of the dimerization domain in the V902L mutant as the critical conformational switch in the cyclase going from the low to high activity state.

Exploring Post-activation Conformational Changes in Pigeon Cryptochrome 4.

A widespread hypothesis ascribes the ability of migratory birds to navigate over large distances to an inclination compass realized by the protein cryptochrome in the birds' retinae. Cryptochromes are activated by blue light, which induces a radical pair state, the spin dynamics of which may become sensitive to earth's weak magnetic fields. The magnetic information is encoded and passed on to downstream processes by structural rearrangements of the protein, the details of which remain vague. We utilize extensive all-atom molecular dynamics simulations to probe the conformational changes of pigeon cryptochrome 4 upon light activation. The structural dynamics are analyzed based on principal component analysis and with the help of distance matrices, which reveal significant changes in selected inter-residue distances. The results are evaluated and discussed with reference to the protein structure and its putative function as a magnetoreceptor. It is suggested that the phosphate-binding loop could act as a gate controlling the access to the flavin adenine dinucleotide cofactor depending on the redox state of the protein.

Introducing Pep McConst-A user-friendly peptide modeler for biophysical applications.

We are introducing Pep McConst-a software that employs a Monte-Carlo algorithm to construct 3D structures of polypeptide chains which could subsequently be studied as stand-alone macromolecules or complement the structure of known proteins. Using an approach to avoid steric clashes, Pep McConst allows to create multiple structures for a predefined primary sequence of amino acids. These structures could then effectively be used for further structural analysis and investigations. The article introduces the algorithm and describes its user-friendly approach that was made possible through the VIKING online platform. Finally, the manuscript provides several highlight examples where Pep McConst was used to predict the structure of the C-terminal of a known protein, generate a missing bit of already crystallized protein structures and simply generate short polypeptide chains.